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THE DIVISIBLE SANDPILE AT CRITICAL DENSITY 


LIONEL LEVINE, MATHAV MURUGAN, YUVAL PERES AND BARIS EVREN UGURCAN 

Abstract. The divisible sandpile starts with i.i.d. random variables (“masses”) 
at the vertices of an infinite, vertex-transitive graph, and redistributes mass by 
a local toppling rule in an attempt to make all masses < 1. The process sta¬ 
bilizes almost surely if m < 1 and it almost surely does not stabilize if m > 1, 
where m is the mean mass per vertex. The main result of this paper is that in 
the critical case m = 1, if the initial masses have finite variance, then the pro¬ 
cess almost surely does not stabilize. To give quantitative estimates on a finite 
graph, we relate the number of topplings to a discrete bi-Laplacian Gaussian 
field. 


1. Introduction 

This paper is concerned with the dichotomy between stabilizing and exploding 
configurations in a model of mass redistribution, the divisible sandpile model. The 
main interest in this model is twofold. First, it is a natural starting place for the 
analogous and more difficult dichotomy in the abelian sandpile model. Second, 
the divisible sandpile itself leads to interesting questions in potential theory. For 
example, under what conditions must a random harmonic function be an almost 
sure constant? (Lemma 5.4 gives some sufficient conditions.) Both the motivation 
for this paper and many of the proof techniques are directly inspired by the work 
of Fey, Meester and Redig [FMR09] . 

By a graph G = (V, E ) we will always mean a connected, locally finite and 
undirected graph with vertex set V and edge set E. We write x ~ y to mean that 
(x, y ) € E. and deg(x) for the number of y such that x ~ y. A divisible sandpile 
configuration on G is a function s : V —> R. We refer to s(x) as an amount of 
‘mass’ present at vertex x; a negative value of s(x) can be imagined as a ‘hole’ 
waiting to be filled by mass. A vertex x € V is called unstable if s(x) > 1. An 
unstable vertex x topples by keeping mass 1 for itself and distributing the excess 
s(x) — 1 equally among its neighbors y ~ x. At each discrete time step, all unstable 
vertices topple simultaneously. (This parallel toppling assumption is mainly for 
simplicity; in Section 2 we will relax it.) The following trivial consequence of 
the toppling rule is worth emphasizing: if for a particular vertex x the inequality 
s(x) > 1 holds at some time, then it holds at all later times. 
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Note that the entire system evolves deterministically once an initial condition 
s is fixed. The initial s can be deterministic or random; below we will see one 
example of each type. Write u n (x) for the total amount of mass emitted before 
time n from x to one of its neighbors. (By the symmetry of the toppling rule, it 
does not matter which neighbor.) This quantity increases with n, so u n | u as 
n f oo for a function u : V —>• [0, oo]. We call this function u the odometer of s. 
Note that if u(x) = oo for some x, then each neighbor of x receives an infinite 
amount of mass from x, so u(y) = oo for all y ~ x. We therefore have the following 
dichotomy: 

Either u(x) < oo for all x € V, 
or u(x) = oo for all x € V. 


In the former case we say that s stabilizes , and in the latter case we say that s 
explodes. 

The following theme repeats itself at several places: The question of whether 
s stabilizes depends not only on s itself but also on the underlying graph. For 
instance, fixing a vertex o, we will see that the divisible sandpile 


s(x) 


1 X J- o 

2 x = o 


stabilizes on G if and only if the simple random walk on G is transient (Lemma 2.9). 

Our main result treats the case of initial masses s(x) that are independent 
and identically distributed (i.i.d.) random variables with finite variance. Write 
Es and Vars for the common mean and variance of the s(x). The mean Es is 
sometimes called the density (in the physical sense of the word, mass per unit 
volume). Because sites topple when their mass exceeds 1, intuition suggests that 
the density should be the main determiner of whether or not s stabilizes: the 
higher the density, the harder it is to stabilize. Indeed, we will see that s stabilizes 
almost surely if Es < 1 (Lemma 4.2) and explodes almost surely if Es > 1 (Lemma 
4.1). Our main result addresses the critical case Es = 1. 


Theorem 1.1. Let s be an i.i.d. divisible sandpile on an infinite, vertex-transitive 
graph, with Es = 1 and 0 < Vars < oo. Then s almost surely does not stabilize. 


Our theme that ‘stabilizability depends on the underlying graph’ repeats again 
in the proof of Theorem 1.1. The proof splits into three cases depending on the 
graph. The cases in increasing order of difficulty are 

• recurrent (Lemma 4.4). Examples: Z, Z 2 . 

• transient with Y^ x ev s(°^ x ) 2 = oo (Section 5.1). Examples: Z 3 ,Z 4 . 

• transient with 9 , (°> x ) 2 < 00 (Section 5.2). Examples: Z d with d > 5. 

Here g denotes Green’s function: g(o, x) is the expected number of visits to x by 
a simple random walk started at o. The reason for the order of difficulty is that 
‘stabilization is harder in lower dimensions,’ in a sense formalized by Theorem 1.2 
below. 
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1.1. Potential theory of real-valued functions. The graph Laplacian A acts 
on functions u : V —>• R by 



( 1 ) 


Using a ‘least action principle’ (Proposition 2.5), the question of whether a divisible 
sandpile stabilizes can be reformulated as a question in potential theory: 

Given a function s : V —>• M, does there exist a nonnegative function 
u : V —>• K such that s + A u < 1 pointwise? 

1.2. Potential theory of integer-valued functions. In the related abelian 
sandpile model, configurations are integer-valued functions s : V —>• Z. We think 
of s(x) as a number of particles present at x. A vertex x € V is unstable if it has 
at least deg(x) particles. An unstable site x topples by sending one particle to 
each of its deg(x) neighbors. This model also has a dichotomy between stabilizing 
(u < oo) and exploding (u = oo), which can be reformulated as follows: 

Given a function s : V —» Z, does there exist a nonnegative function 
u : V —>• Z such that s + A u < deg —1 pointwise? 

The restriction that u must be integer-valued introduces new difficulties that 
are not present in the divisible sandpile model. The first step in the proof of 
Theorem 1.1 is to argue that if Es = 1 and s stabilizes, then it necessarily stabi¬ 
lizes to the all 1 configuration. This step fails for the abelian sandpile except in 
dimension 1. Indeed, a result analogous to Theorem 1.1 does hold for the abelian 
sandpile when the underlying graph is Z [FMR09, Theorem 3.2], but no such result 
can hold in higher dimensions: The density Es alone is not enough to determine 
whether an abelian sandpile s on Z rf stabilizes, if d < Es < 2d — 1 (see [FR05, 
Section 5], [FMR09, Theorem 3.1] and [FLP10, Proposition 1.4]; the essential idea 
in these arguments arose first in bootstrap percolation [Ent87, Sch92]). 

We would like to highlight an open problem: Given a probability distribution p 
on Z (say, supported on {0,1,2,3,4} with rational probabilities) is it algorithmi¬ 
cally decidable whether the i.i.d. abelian sandpile on Z 2 with marginal p stabilizes 
almost surely? 

1.3. Quantitative estimates and bi-Laplacian field. For a finite connected 
graph G = (V, E), the divisible sandpile s : V —>• M stabilizes if and only if 
YIxgV s ( x ) — M- Our next result gives the order of the odometer in a critical 
case when this sum is exactly |V|. Specifically, to formalize the idea that ‘stabi¬ 
lization is harder in lower dimensions,’ we take an identically distributed Gaussian 
initial condition on the discrete torus Z^, conditioned to have total mass n d . The 
expected odometer can be taken as an indication of difficulty to stabilize: How 
much mass must each site emit on average? According to equation (2) below, the 
expected odometer tends to oo with n in all dimensions (reflecting the failure to 
stabilize on the infinite lattice Z d ) but it decreases with dimension. 
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Theorem 1.2. Let (cr(x)) xgZ d be i.i.d. iV(0,1), and consider the divisible sandpile 

s dtn (x) = 1 + cr(x) -—^2 cr(y). 

U y&K 

Then s d n : —>• R stabilizes to the all 1 configuration, and there exists a constant 

Cd such that the odometer u dyTl satisfies 


c d 1 fid(n) < E u dt n(x) < C d (j) d (n 

for all n > 2, where cf> d is defined by 



"n 3 / 2 , 

d = 1 


n, 

d = 2 

fidfi) := < 

n 1 / 2 , 

d = 3 


log n, 

II 


(log n) 1 / 2 , 

d > 5. 


The first step in computing these orders is proving an equality in law between 
the odometer u diU and a certain ‘discrete bi-Laplacian Gaussian field’ shifted to 
have minimum value 0. This equality in law actually holds for any finite connected 
graph, as detailed in the next proposition. 

Proposition 1.3. Let G = (V,E) be a finite connected graph. Let (cr(x)) xe y be 
i.i.d. N(0,1), and consider the divisible sandpile 

s(x) = i + cr{x) - a (y)- 

11 yev 

Then s stabilizes to the all 1 configuration, and the distribution of its odometer 
u : V —>• [0, oo) is 

(u(x)) xe v = (r/(x) - min rj) xeV 

where the rj{x) are jointly Gaussian with mean zero and covariance 

= degWdegfa) Y.sa, X )a{z,V) 

where g is defined by g(x,y) = YlzeV 9 z i x -> V) an d iff- V) the expected num¬ 
ber of visits to y by the simple random walk started at x before hitting z. 

Proposition 1.3 suggests the possibility of a central limit theorem for the di¬ 
visible sandpile odometer on Z^: We believe that if cr is identically distributed 
with zero mean and finite variance, then the odometer, after a suitable shift and 
rescaling, converges weakly as n —» oo to the bi-Laplacian Gaussian field on M d . 
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1.4. Proof ideas. By conservation of density (Proposition 3.1) the assumption 
Es = 1 implies that if s stabilizes then it must stabilize to the all 1 configuration, 
so that the odometer u satisfies 

A u = 1 — s. (3) 

where A is the Laplacian (1). This relation leads to a contradiction in one of three 
ways: 

• If G is recurrent (examples: Z, Z 2 ), then 1 +f55 0 does not stabilize (Lemma 2. 
By resampling the random variable s(o) we derive a contradiction from (3). 

• If G is simply transient (examples: Z 3 ,Z 4 ) then we can attempt to solve 
(3) for u, writing 

u (v) = ^9(x,y)(s{x) - 1) (4) 

X 

where g is Green’s function. The sum on the right side diverges a.s. if taken 
over all x € V (since g(-,y) is not square-summable) but we can stabilize 
s in nested finite subsets V n ^ V instead. The corresponding finite sums, 
suitably normalized, tend in distribution to a mean zero Gaussian by the 
Lindeberg central limit theorem, contradicting the nonnegativity of u(y). 

• If G is doubly transient (example Z rf for d > 5) then the right side of (4) 
converges a.s.. The difference between the left and right sides is then a 
random harmonic function with automorphism-invariant law. The proof 
is completed by showing that under mild moment assumptions any such 
function is an almost sure constant (Lemma 5.4). 

1.5. Related work. The divisible sandpile was introduced in [LP09, LP10] to 
study the scaling limits of two growth models, rotor aggregation and internal 
DLA. The divisible sandpile has also been used as a device for proving an exact 
mean value property for discrete harmonic functions [JLS13, Lemma 2.2], These 
works focused on sandpiles with finite total mass on an infinite graph, in which 
case exploding is not a possibility. In the present paper we expand the focus to 
sandpiles with infinite total mass. 

The abelian sandpile has a much longer history: it arose in statistical physics 
as a model of ‘self-organized criticality’ (SOC) [BTW87, Dha90]. The dichotomy 
between stabilizing and exploding configurations arose in the course of a debate 
about whether SOC does or does not involve tuning a parameter to a critical value 
[FR05, MQ05]. Without reopening that particular debate, we view the stabiliz¬ 
ing/exploding dichotomy as a topic with its own intrinsic mathematical interest. 
An example of its importance can be seen in the partial differential equation for 
the scaling limit of the abelian sandpile on Z 2 , which relies on a classification of 
certain ‘quadratic’ sandpiles according to whether they are stabilizing or exploding 
[LPS12], 

The Gaussian vector rj in Proposition 1.3 can be interpreted as a discrete bi- 
Laplacian field. In Z d for dimensions d > 5, Sun and Wu construct another discrete 
model for the bi-Laplacian field by assigning random signs to each component of 
the uniform spanning forest [SW13]. 
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2. Toppling procedures and stabilization 

In this section G = (V, E) is a locally finite, connected, undirected graph. 
Denote by X = M' the set of divisible sandpile configurations on G. 

Definition 2.1. Let T C [0,oo) be a well-ordered set of toppling times such that 
0 € T and T is a closed subset of [0, oo). A toppling procedure is a function 

T x V —» [0, oo) 

(t,x) ufix) 

satisfying for all x € V 

(1) u 0 (x) = 0. 

(2) ut 1 (x) < ut 2 {x ) for all t\ < t 2 - 

(3) If t n f t, then u tn (x) | u t (x). 

In the more general toppling procedures considered by Fey, Meester and Redig 
[FMR09], the assumption that T is well-ordered becomes a “no infinite backward 
chain” condition, but we will not need that level of generality. See Examples 1-3 
below for the three specific toppling procedures we will use. 

The interpretation of a toppling procedure is that starting from an initial con¬ 
figuration s € X, the total mass emitted by a site x € V to each of its neighbors 
during the time interval [0,t] is Ut(x), so that the resulting configuration at time 
t is 

s t = s + A ut 

where A is the graph Laplacian (1). 

For a € M write a + = max(a, 0). For t G T write t~ := sup{r € T : r < t}. 
Note that t~ € T since T is closed. 


Definition 2.2. A toppling proceedure u is called legal for initial configuration s 
if 


ut(x) - u t -(x) < 


(s t -(x) - 1)+ 
deg(x) 


for all and all t G T \ {0}. 

Thus, in a legal toppling procedure, a site with mass < 1 cannot emit any mass, 
while a site with mass > 1 must keep at least mass 1 for itself. 

Definition 2.3. A toppling procedure u is called finite if for all x € V we have 

Uoo(x) := lim uAx) < oo 

t->supT 


and infinite otherwise. The limit exists in [0, oo] since ut{x) is nondecreasing in t. 
Note that if u is a finite toppling procedure, then the limit 
s oo := lim st = s + lim A ut 

t—y sup T t—>-sup T 


exists and equals s + Au^. 
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Definition 2.4. Let s E X. A toppling procedure u is called stabilizing for s if u 
is finite and Soo < 1 pointwise. We say that s stabilizes if there exists a stabilizing 
toppling procedure for s. 

Throughout this paper, all inequalities between functions hold pointwise, and 
we will usually omit the word “pointwise.” 

A basic question arises: For which s E X does there exist a stabilizing top¬ 
pling procedure? For instance, one might expect (correctly) that there is no such 
procedure for s = 2. We can rephrase this question in terms of the set of functions 

F s := {/ . v -» R | / > 0 and s + A/ < 1}. 

If u is a stabilizing toppling procedure for s, then Uoo E F s . Conversely, any / E F s 
arises from a stabilizing toppling procedure for s simply by setting T = {0,1} and 
u\ = /. Therefore s stabilizes if and only if F s is nonempty. 

Proposition 2.5. (Least action principle and abelian property) Let and 

let l be a legal toppling procedure for s. 

(i) For all f € F s , 

loo </. 

(ii) If u is any stabilizing toppling procedure for s, then 

IoG if • 

(iii) If u is any legal stabilizing toppling procedure for s, then for all x E V, 

Uoo(x) = inf{/(.x) | / € F s } (5) 

In particular, u a0 and the final configuration 

Sqq — s T Auqo 

do not depend on the choice of legal stabilizing toppling procedure u. 

Proof, (i) For y E V let r y = inf{f E T : itiy) > /(y)}, and suppose for a 
contradiction that r y < oo for some y E V. Since T is well-ordered, the infimum 
is attained. Moreover, r y := sup{t E T : t < r y } < r y by assumption (3) of 
Definition 2.1. 

Let r = infygy r y . Since T is well-ordered, r = r y for some y E V. Now at time 
t~ , since u is legal for s, 

S T - (y)> 1 + deg(y) (4(y) - 4- (y)) > 1 + deg (y)w(y) 
where w = f — I T ~ . On the other hand, 

s T -(y) = (s +A e T -)(y) 

= (s + Af)(y) - A w(y) 

< 1 + deg (y)w(y) - E w(x). 

x^y 

It follows that Yhx~ y w (x) < 0. But for all x E V we have t~ < t < t x , so 
w(x) > 0, which yields the required contradiction. 

Part (ii) follows from (i), using / = u 0Q . 
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Part (iii) also follows from (i): 
infimum for all x € V. 


the function Uoo simultaneously attains the 

□ 


Whenever we need to fix a particular toppling procedure, we will choose one of 
the following three. 


Example 1. (Toppling in parallel) Let T = N. At each time t € N, all unstable 
sites of St- 1 topple all of their excess mass simultaneously: For all x € V, 


u t {x) - u t -\{x) 


(gt-i(a) - 1) + 
deg(x) 


This u is a legal toppling procedure. Two further observations about u will be 
useful in the proof of Lemma 2.7 below. First, if s stabilizes, then u is finite by 
Proposition 2.5(h). Second, whenever u is finite, u is stabilizing for s: indeed, if 
Uoo(x) = J2tm( s t( x ) - X ) + < °°> we have ( s t( x ) - 1) + ->■ o as t -A oo and 

hence Soo(x) < 1. 


Example 2. (Toppling in nested volumes) Let V\ C Vz C ... be finite sets with 
Un>i Ki = V . Between times n — 1 and n we topple in parallel to stabilize all 
sites in V n : Formally, we take T to be the set of all rationals of the form n — £ for 
positive integers n and k. For n > 1 and k > 1 we set 


and 


u n i{x)~u n 1 (x) 

fc k — 1 


(s n _ l(x) - 1 )+ • 1 X £V„ 

k — 1 

deg(x) 


u n (x) = lim u i(x). 

k ->00 n k 


Example 3. (Toppling in two stages) In this procedure we are given a decomposi¬ 
tion of the initial configuration into two pieces 


where s 1 stabilizes and s 2 > 0. In the first stage we ignore the extra mass s 2 and 
stabilize the s 1 piece by toppling in parallel at times 1 — ^ for positive integers k , 
obtaining 

si = s + = s^ + s 2 . 

The condition that s 1 stabilizes ensures u ^ < 00 , and the condition s 2 > 0 ensures 
that all topplings that are legal for s 1 are also legal for s. Now we topple s^ + s 2 
in parallel at times 2 — ^ for positive integers k. 

Now we come to a central definition of this paper. Let s € X. 


Definition 2.6. The odometer of s is the function : V —>• [0, 00 ] of (5). If s 
stabilizes, then its stabilization is the configuration 

Sqo — s T Ai/qq. 
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If s stabilizes, then its odometer Uoo(x) is the total amount of mass sent from x 
to one of its neighbors, in any legal stabilizing toppling procedure for s. If s does 
not stabilize, then = oo: The odometer is defined as a pointwise infimum, with 
the usual convention that the infimum of the empty set is oo. Next we observe 
that the odometer can also be expressed as a pointwise supremum. 

Lemma 2.7. Let Uqo be the odometer of s € X. Then for all x E V, 

u 00 (x) = sup{£oo(x) | £ is a legal toppling procedure for s}. (6) 

Proof. Denote the right side of (6) by L(x). By Proposition 2.5 (i), L < Uoa. To 
prove the reverse inequality we will use a particular legal l, the parallel toppling 
procedure of Example 1. There are two cases: First, if this l is finite, then i is 
stabilizing as well as legal, so L > = Uoo by Proposition 2.5(iii). 

Second, if £ is not finite, then £ 00 [o) = oo for some o £ V. Then for any neighbor 
x ~ o, we have 

£t+i{x) > £t(o) + s(x) - 1 

and the right side tends to oo with t, so £oo{x) = oo. Since the graph G is 
connected it follows that £oo = oo. In this case, both L and Uqo are identically 
oo. □ 


We pause to record several equivalent conditions for s stabilizing. 


Corollary 2.8. Let s € X have odometer Uoq. The following are equivalent. 

(1) There exists a legal stabilizing toppling procedure for s. 

(2) There exists a stabilizing toppling procedure for s. 

(3) T s + §. 

(4) Uoo(x) < oo for all x € V. 

(5) Every legal toppling procedure for s is finite. 

(6) The parallel toppling procedure for s is finite. 

Proof. The implications (1) (2) (3) => (4) and (5) => (6) are immediate from 

the definitions. The implication (4) => (5) follows from Lemma 2.7. Finally, let 
u be the parallel toppling procedure for s. If u is finite, then u is both legal and 
stabilizing for s, which shows the remaining implication (6) => (1). □ 


Denote by P Q the law of the discrete time simple random walk (Xj)j e ^ on 
G started at Xq = o. Writing pj(x) = JP 0 (W/' = x)/deg(x), observe that the 
Laplacian (1) of pj satisfies 


Xpj{x) 


E 

yr^X 


f FpjXj = y) 
V deg (y) 


F 0 (Xj = x) \ 
deg(x) J 


= ^2 p o(Xj = y, X j+ 1 = x) - P 0 (Xj = x) 

yr^X 

= deg(x)(p j+1 (x)-pj(x)). (7) 


The next lemma will play an important role in the proof of the recurrent case of 
Theorem 1.1. It relates the stabilizability of a particular configuration s to the 
transience of the simple random walk. 
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Lemma 2.9. Fix o £ V and (3 > 0. The divisible sandpile 


s(x) 


1 , X 7 ^ o 

1 + (3, x = o 


stabilizes on G if and only if the simple random walk on G is transient. 


Proof. We compute the parallel toppling procedure ut and the configuration st at 
time t of Example 1. Setting pt(x) = ¥ a (X t = x)/ deg(x), let us show by induction 
on t that 

t -1 

M x ) = P^2pji x ) ( 8 ) 

j =0 

and 

s t (x) = 1 + f3 deg (x)pt(x) (9) 

for all x £ V and all t £ N. Indeed, if (8) holds at time t, then summing (7) we 
obtain 

st = s + A ut = s + /3 deg (x)(p t - po) = 1 + (3 deg (x)p t 
so (9) holds at time t, whence 

u t+ i(x) - u t (x ) = ^ = /3p t (x) 

deg(x) 

so (8) holds at time t + 1, completing the inductive step. 

Taking 1 1 oo in (8) yields Uoo(x) = (3g(o,x)/ deg(x) where 

OO 

9(o, x) = ^ P 0 (Xj = x) (10) 

3 =0 

is the Green function of G, which is finite if and only if G is transient. 

If G is transient, then the parallel toppling procedure u is finite and s+Ai/qo = 1, 
so s stabilizes. If G is recurrent, then the parallel toppling procedure is infinite, 
so s does not stabilize by Corollary 2.8. □ 


3. Conservation of density 

In this section we assume that G = ( V, E ) is vertex-transitive, and we fix 
a subgroup T of Aut(G') that acts transitively: for any x,y £ V there is an 
automorphism a £ T such that ax = y. For the rest of the paper, we assume 
o £ V be an arbitrary fixed vertex. Write X = M' (viewed as a measurable space 
with the Borel u-field) and T a : X —>• X for the shift ( T a f)(x ) = f(a~ 1 x). Let 
P be a probability measure on X satisfying E|s(o)| < 00 . We assume that P is 
T-invariant; that is, if s has distribution P then T a s has distribution P for all 
a £ T. 

Proposition 3.1. (Conservation of Density) If ¥ is T-invariant and F{s stabilizes} 
1, then the stabilization s ^ satisfies 

Psoo(o) = Es(o). 
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Fey, Meester and Redig [FMR09] used an ergodic theory argument to prove the 
conservation of density for the abelian sandpile on when T is the group of trans¬ 
lations. They considered only nonnegative initial conditions: s > 0. We will give 
an elementary proof which starts in the same way, by choosing an automorphism- 
invariant toppling procedure (toppling in parallel) and uses the following observa¬ 
tion about averages of uniformly integrable random variables. 

Lemma 3.2. If X is a random variable with E|X| < oo, and X\, X 2 , ■ ■ ■ is any 
sequence of random variables such that X % = X for all i, then the family 

S = {a\Xi + ... + a k X k | fc > 1, cti > 0, Qj = 1} 
is uniformly integrable. 

Proof. We may assume X > 0. Given e > 0, we must show that there is a <5 > 0 
so that for any set A with P(A) <6 we have ECKl^) < e for all Y e S. 

Since EX < 00 , there is such 6 for X itself: choose M so that E(X1 {x >M }) < e 
and then set 5 = P(X > M ). Now for any set A with P(A) < 5 and any Y = 

y. diXi € s 

k 

E[Y1 A ] =J2 a i EPQIa] < e 

2—1 

so the same 5 works for all Y £ S. □ 

Proof of Proposition 3.1. We topple in parallel: at each time step t = 0,1,..., 
each site x € V distributes all of its excess mass = ( st(x ) — 1) + equally 

among its r neighbors where r is the common degree of all vertices in G. The 
resulting configuration after t time steps is 

st = s 0 + A u t . 

where ut = r~ 1 (cro + ... + at- 1 ). Since P{s stabilizes} = 1 we have st(o ) —>• 
Sco(o), a.s.. We will show that the random variables at(o) for t £ N are uniformly 
integrable. To finish the proof from there, note that the law of ut is T-invariant, 
so E utfx) = E ut(y) for all x,y € V. In particular, EA ut.(o) = 0 and hence 
Esi(o) = Eso(o) for all t < 00 . Since st > min(so, 1) the uniform integrability of 
at(o) implies that of st(o), so we conclude Esoo(o) = Eso(o). 

At time step t the origin retains mass < 1 and receives mass at~i{y)/r from 
each neighbor y, so 

»,(»)< i + X^ 

x~o 

hence 

at(o) < -^at-ijx). 

x~o 

Inducting on t we find that 

a t {o) < Y t := ^ a t (x)a 0 {x) 
xev 



12 


LEVINE, MURUGAN, PERES, UGURCAN 


where o*(x) is the probability that a t-step simple random walk started at o £ V 
ends at x. By Lemma 3.2 the random variables Yt are uniformly integrable, which 
completes the proof. □ 

We remark that the above proof also applies to the abelian sandpile to show 
[FMR09, Lemma 2.10], by taking &t(x) = \_st(x) + /r\, the number of times x 
topples at time t. 

4. Behavior at critical density 

In this section G = ( V , E) is an infinite vertex-transitive graph, o £ V is a fixed 
vertex, T C Aut(G) is a group of automorphisms that acts transitively V, and P is 
a T-invariant and ergodic probability measure on X = R' with E|s(o)| < oo. The 
event that s stabilizes is T-invariant, so it has probability 0 or 1 by ergodicity. 

Lemma 4.1. //Es(o) > 1, then P{s stabilizes} = 0. 

Proof. If s stabilizes, then by conservation of density (Proposition 3.1), we have 
E Soo (o) = Es(o). Since the final configuration s 00 is stable we have Soo(o) < 1, so 
Es(o) <1. □ 

Lemma 4.2. //Es(o) < 1 then P{s stabilizes} = 1. 

Proof. We will show the contrapositive. We will use the observation made in the 
introduction that if St(o) > 1 for some time t then st(o) > 1 for all T > t. 

Consider first the case that s is bounded below: P(s(o) > M ) = 1 for some 
M £ (—oo,0]. Define ut and st by toppling in parallel as in §3. Supposing that 
P{s stabilizes} = 0, we have ut(o ) > 0 for some sufficiently large t, a.s.; this fact is 
contained in the proof of Lemma 2.7. Since ut(o) > 0 implies st(o) > 1 and hence 
st(o ) > 1 for all T > t, we have 

P jliminf St(o) > l| = 1. 

Since st > min(s, 1) > M, by Fatou’s lemma E (liminf sj(o)) < Es(o), which 

V t—>oo J 

shows Es(o) > 1 as desired. 

Now we use a truncation argument to reduce the general case to case where s is 
bounded below. Choose sufficiently small M £ (—oo,0] such that Es(o)1<,( 0 )xm < 
1. By the previous case the configuration s1 s >m stabilizes almost surely. Since 
s < s1 s >m, we have that s stabilizes almost surely. □ 

If s stabilizes then the odometer of s is the function u 0Q : V —» [0,oo), where 
u is any legal stabilizing toppling procedure for s. 

Lemma 4.3. Ifs has T-invariant law P and¥{s stabilizes} = 1, then the odometer 
Uoo has r -invariant law. Moreover, ifKs(o) = 1 then s 00 = 1 and Au^ = 1 — s. 

Proof. We use (5) along with the fact that A commutes with T a : if s + A/ < 1 
then T a s + A (T a f) = T a (s + A/) < 1, so if u 0 0 is the odometer for s then T a Uoo 
is the odometer for T a s. 

By conservation of density (Proposition 3.1), Esoo(o) = 1. Since Sqo < 1 it 
follows that Sqo = 1, and hence Auoo = 1 — s. □ 
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In the case that G is recurrent we can prove our main theorem with no moment 
assumption, using an “extra head” construction. 

Lemma 4.4. Suppose that G is an infinite, recurrent, vertex-transitive graph and 
{s(x)} Ig y are i.i.d. with Es(x) = 1 and P{s(x) = 1} ^ 1. Then 

P{s stabilizes} = 0. 

Proof. Let S C X denote the set of configurations that stabilize. By ergodicity, 
P(S) € {0,1}. We will show that if P(S') = 1, then the graph G must be transient. 
Let s : V —> R denote an i.i.d. configuration with the given distribution such that 
Es = 1, and fix a vertex o E V. We create a new i.i.d. configuration s' with the 
same law as s, by independently resampling s(o) from the same distribution. Then 
s' = s + /35 0 , where 5 0 (x) = \{x = o} and fi is a mean zero random variable. Since 
Vars > 0, we have that P(/3 > 0) > 0. Using P(S') = 1 along with Lemma 4.3, 
there exist s : V —>• R and /3 > 0 such that s and s + /35 0 both stabilize to the all 
1 configuration. By toppling s + f35 0 in two stages (Example 3), it follows that 
1 + f35 0 stabilizes to 1, so G is transient by Lemma 2.9. □ 

In the preceding lemma the hypothesis that s is i.i.d. can be substantially weak¬ 
ened: The proof uses only the fact that with positive probability, the conditional 
distribution of s(o) given {s(x)} XJ t 0 is not a single atom. 

5. Proof of Theorem 1.1 

5.1. Singly transient case. Recall Green’s function (10). In this section we 
assume that g(o, y) < oo for all y € V but 

X^(o,y) 2 = oo. (11) 

yev 

Define V n = {x E V : d{x, o) < n} where d is the graph distance. Then Vi C 
Vi C ... are finite sets with U n >i Hi = H. Let g n {x,y) be the expected number 
of visits to y by simple random walk started at x and killed on exiting V n . By the 
monotone convergence theorem, for fixed x, y € V we have 

9n(x,y) t g{x,y) 

as n —>• oo. In particular, setting 

( V 2 

Vn := 9n{o,xf 

\y&v n ) 

we have u n t oo as n —> oo by (11). 

The proof of the following lemma is inspired by [FR05, Theorem 3.1] and 
[FMR09, Theorem 3.5]. 

Lemma 5.1. Let G = {V. E) be singly transient (11), 

{s(x)}a; e y be i.i.d. with Es = 1 and Vars < oo. If 

— 9 n(o,x)(s(x) - 1) 

^n 


vertex transitive graph. Let 


14 


LEVINE, MURUGAN, PERES, UGURCAN 


converges in distribution as n —>• oo to a nondegenerate normal random variable Z, 
then P{s stabilizes} = 0. 


Proof. Assume to the contrary that s stabilizes a.s.. Let u be the nested volume 
toppling procedure (Example 2). Then for each n € N, 

s + A u n = on V n (12) 

with i n < 1 on V n . Equation (12) can be rewritten as 

Un(y) = Y 9 n{x,y){s(x ) - 1) +r _1 Y 9n(x,y)( 1 -in(x)) 

XGVn x£Vn 


where r is the common degree (both sides have Laplacian f n — s in V n and vanish on 
V^). Observe that the second term is a nonnegative random variable. Therefore, 
for any e > 0 


F{u n (o) > eis n } > P 



Y 9 n(o,x)(s(x) - 1) 
x£V n 


> re 


(13) 


Since u is a legal toppling procedure and we have assumed that s stabilizes a.s., 
we have rtoo(o) < oo, a.s.. Now since u n (o) t Uoo(o) and v n t oo, the left side of 
(13) tends to zero as n —» oo. However, the right side tends to a positive limit 
P(Z > re) > 0, which gives the desired contradiction. □ 


To complete the proof of Theorem 1.1 in the singly transient case, it remains 
to show that x Yxev 9 n(o, x)(s(x) — 1) converges in distribution to a nondegen¬ 
erate normal random variable. We show this using Lindeberg-Feller central limit 
theorem as follows. 


Lemma 5.2. Let {X n \ : n > 1; i = 1,..., k n } be a triangular array of identically 
distributed random variables such that for each n € N, {X m : i = 1,..., k n } is 
independent. Assume that E[X n] = 0 and E[X\ X ] = 1. Let a n k > 0 be such that 
YkL i a nk = 1 an d li m n->oo^n = 0 where b n = maxi<fc<fc n a n k- Then the sequence 
Y n = YY a nkX n k converges in distribution to standard normal random variable 
as n —>• oo. 


Proof. By Lindeberg-Feller central limit theorem [Bil95, Theorem 27.2] it suffices 
to check the Lindeberg condition: 

kn 

|>e] = 0 

k =1 

for all e > 0. Since E[-X£ fc l| 0nfcX „ h |> e ] < E [ x n 1 \x 11 \>e/bJ and Yk a lk = 1 . we have 

kn 

} i ^Jl E i a 2 nk X nk 1 KkX nk \>e\ < n 1 ™ 0 E [ X n 1 ^i 1 |> 6 /b„] = 0 
k =1 


because liiUn^oo b n = 0 and = 1. 


□ 
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Lemma 5.3. Let G = (V. E) be a singly transient (11), vertex transitive graph. 
Suppose {s(x)}a; e v are i.i.d. random variables with Es = 1 and 0 < Vars < oo. 
Then P{s stabilizes} = 0. 


Proof. Since g(o , o) < oo and v n t oo, we have 

r g n (o,x ) g n {o,o ) g(o,o ) 

lim max-= lim -< lim -= (J. 

n—>-oo xeVn 


For n € N and x € V n , define X nx = (s(x) — l)/\/Var s and a 
b n = ma x x£ y n a nx . By Lemma 5.2, we have that 

— Y 9 n{o,x)(s(x) - 1) 

Vn .. 

XGVn 


n—» oo 


ffn(o,x) 

Vn 


and 


converges in distribution to a normal random variable with mean 0 and variance 
Var(s). Lemma 5.1 implies the desired conclusion. □ 


5.2. Doubly transient case. We start by outlining our proof strategy. Recall 
from Lemma 4.3 that if s stabilizes with Es(o) = 1 then it must stabilize to the 
constant configuration s 00 = 1. In particular, the odometer u 00 satisfies Artoo = 
1 — s. In Lemma 5.5 below, by convolving s — 1 with Green’s function we can 
build an explicit function v with Laplacian 1 — s; the convolution is defined almost 
surely provided that 

X^5(o>y) 2 < oo. (14) 

y&V 

(This condition says that the expected number of collisions of two independent 
random walks started at o is finite. The essential feature of Green’s function here 
is of course that 

Aff(-,y) = -rS y (15) 

where r is the common degree of all vertices of G, and 5 y (x) = l{x = y}.) 

Having built the function v, the difference v — Uoo is then a random harmonic 
function with T-invariant law, which must be an almost sure constant by the 
following lemma. 

Lemma 5.4. (Harmonic Functions With Invariant Law) Let T be a group of 
automorphisms of G that acts transitively on the vertex set V. Suppose that h : 
V —>• M has T-invariant law and Ah = 0. If h can be expressed as a difference of 
two functions h = v — u where u > 0 and sup^gy Ev(x) + < oo, then h is almost 
surely constant. 

Proof. Let (X n ) n > 0 be simple random walk on G started at Xo = o. Although h is 
harmonic, h(X n ) need not be a martingale since it need not have finite expectation. 
But for any a € R, since h = v — u and u > 0, the truncation 

M n := a + {h(X n ) - a) + 

<a + {v(X n ) - a) + 
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has finite expectation. Since the function t i-A (t — a) + is convex, x i-A (h(x) — a) + 
is subharmonic, so M n is a submartingale: 

E[M n+1 \h, X 0 , ...,X n ] = - r £ (a + (>*H - «) + ) > M n■ 

W~X n 

A submartingale bounded in L 1 converges almost surely [Wil91, 11.5]. Since this 
holds for any a € M it follows that h(X n ) converges almost surely. But since h has 
T-invariant law, h(X n ) is a stationary sequence, so the only way it can converge 
a.s. is if h(X n ) = h(X o) for all n. Since G is connected, for any vertex x E V there 
exists n with P(X n = x) > 0, so h(x) = h(o). □ 


Lemma 5.5. Let T be a countable subgroup of Aut(G) which acts transitively on 
the vertices of V. Suppose {a(x)} x£ v i.i.d. random variables with Ecr(x) = 0 
and 0 < Varcr(O) < oo Let 2 / 1 , 2 / 2 , - - - be an enumeration of the vertex set of G. 
For a£f, let 


V a (x) := ~^2g{x,ayi)a(ayi). 

2—1 

If G is doubly transient (14), then the following hold almost surely. 

(a) The series defining v a (x) converges for all x € V,a: € T. 

(b) v a = vy for all a, /3 € T. 

(c) Av e = —a where e denotes the identity automorphism. 

(d) v e has r -invariant law, that is T a v e = v e for all a € T. 

(e) v e is unbounded above and below. 


(16) 


Proof, (a) Each term in the series (16) is independent, and the i-tli term has 
variance g(x,ayi) 2 U where U = r~ 2 Eer(0) 2 . By (14) the sum of these variances is 
finite, which implies that the series converges a.s.. 

(b) Note that there is something to check here because the series defining 
v a is only conditionally convergent. Fix x € V and denote the partial sum 
r -1 X^r=i 9 , ( x ) a yi) a i a yi) by v a , n . We compare v a)n and v e>n where e is the identity 
element of T. Given e > 0, choose n\ such that 

OO 

U y{x,yi ) 2 < e 3 . 

i=ni 

There exists N\ depending on m and a such that 

{yi : i = 1, 2,... , m} C {ay t :i = 1,2,..., Ah}- 

This implies that for any N > N\ we have Var (v a ,N - v e>N ) < 2e 3 . By Cheby- 
shev’s inequality, 

P (\v a ,N - v e ,N | > e) < e" 2 Var ( v a ,N - v e ,N) < 2e 

for all N > N\. By part (a) there a.s. exists A^ > Ni such that m.cQt(\v a)n — 
Uq,(.t)|, \v etn — v e (x)\) < e/3 for all n > AT 2 . By the triangle inequality it follows 
that 

P (V Q (x) - v e (x)\ > |) < 2e. 
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Since e was arbitrary we obtain v a (x) = v e (x), a.s.. By taking countable intersec¬ 
tions we have that v a = v e , a.s.. 

(c) Using (15), we compute 


Av e (x) 


~ -v e (x)) 

W~X 

j / O v 

EE g{w,yi)a{yi) -^ j g{x,y i )a{y i ) 

w~x \ i=l i =1 

\ ^ 

(. g(w,yi) - g{x,yi)) 


i= 1 
oo 


5>(lfc)l{w = 


2=1 


-o-(x). 


(d) To show that v e has T-invariant law, write v e = to make the dependence 
on the initial configuration a explicit. We have for all a € T 

OO 

T a v°(x) = 22g(a~ 1 x 1 yi)a{y i ) 

2=1 

OO 

= ^2g(x,ayi)a(yi) 

2=1 

OO 

= ^2g{x,ayi)(T a cr)(ayi) 

2=1 

= Va aa {x) 

= vj aa (x) 


where in the last equality we have used part (b). Hence T a v a e = = v a e since 

(7 has T-invariant law. 

(e) To show that v e is almost surely unbounded below, we use the assumption 
that <r(o) has zero mean and positive variance, which implies that 

P(cr(o) < —6) > p 

for some p, 5 > 0. Since g(o , y) = oo and 9(°i v) 2 < °°) we can choose 

N large enough so that 

N oo 

(5r_1 ^2 9(o, yi) > 2 M, u "22 g{o,yi) 2 < 1 . 

i=l i=JV+l 

By Chebyshev’s inequality 

X j 

P ( r_1 5Z g(x,yi)a{ yi ) > M) < —j. 

i=N+l 
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On the event that a (y*) < —<5 for all 1 < i < N and r 1 YaLn+i M 

we have v e {o) < —2 M + M. Since the random variables a{yi) are i.i.d., we obtain 

F(v e (o)<-M)>p N ( l-^)>0. 

Since v e is stationary, we have P(inf v e < — M) € {0,1} by ergodicity. Since this 
probability is > 0, it must be 1. Since M was arbitrary, v e is a.s. unbounded 
below. 

A similar argument shows that v e is also a.s. unbounded above. □ 

Lemma 5.6. Let G = (V,E) be a doubly transient (14), vertex transitive graph 
and let {s(x)}a; e y be i.i.d. random variables with Es = 1 and 0 < Vars < oo. 
Then P{s stabilizes} = 0. 

Proof. Assume to the contrary that s stabilizes a.s. with odometer Uoq. Let T be 
a countable subgroup of Aut(G') that acts transitively on the vertices of G. Let 
v : V —>• M be given by equation (16) with a = e and a = s — 1. By Lemma 5.5, v 
has T-invariant law and An = 1 — s. Therefore 

h := v - Uoo 

has T-invariant law and h is harmonic: Ah = 0 on V. Further, by Fatou’s lemma, 
Eu(x) 2 < deg(x) -2 Var(s) g(x, yf) 1 < oo. Lemma 5.4 now implies that h is 
almost surely constant. This contradicts Lemma 5.5 because rioo > 0 and v is 
almost surely unbounded below. □ 

6. Stabilizability of Cones 

Until now we have mainly been concerned with the stabilizability of random 
initial configurations. In this section we examine stabilizability of a few determin¬ 
istic configurations on the square grid Z 2 . We present two examples, one of which 
stabilizes. 

Lemma 6.1. Define G\ = {(x,y) € Z 2 : x > 0, \y\ < a:}. Then the configuration 
so = (1 T ajlc’i does not stabilize for a > 0. 

Proof. By least action principle (Proposition 2.5), it suffices to show the existence 
of an infinite legal toppling procedure. 

Let Uk denote the parallel toppling procedure of Example 1 where k € N. Let 
C = {(x,y) £Z 2 : x > 0, \y\ < x}. Let (X n ,Y n ) denote the simple random walk 
on Z 2 and let N denote the stopping time N = min{n > 0 : X n = |W,|}. As in the 
proof of Lemma 2.9, we keep track of the mass from each (x, y) E C to obtain 

k 

M 1,0) ^(x,y)(( X N, Y N) = (0,0), N = l ). 

( x,y)GC 1=0 

for all k G N. To see this note that N is the exit time of the set C and the only 
way to exit C at (0,0) is from (1,0). As a result, we have 

Uoc(l,0)>j P( X iV) 

(x,y)€C 
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where p(x, y) 
given by 


p(-£, 2 /) = (0)0))- Consider the function q : Z 2 —>• [0,1] 


q(x,y) 


p(\x\,y) if |y| < M, 

< -p(\y\,x) if |x| < |y|, 
0 otherwise. 


Note that A q = <5( 0 ,i) + <5(o,-i) — ^(i,o) — ^(-i,o)- Let g(x,y ) denote the potential 
kernel in Z 2 defined by 


g(x,y) = X][ p (oh)( x « = ( x ,y)) - F (o,o)pfn = (o, o))] 

n=0 


where (A n ) ng f*j denotes the simple random walk on Z 2 . Although the simple 
random walk on Z 2 is transient, it turns out that the sum defining g is absolutely 
convergent. By standard estimates on g (See [Law96, Chapter 1]), we know that 
g has sub-linear (logarithmic) growth. This combined with [HS93, Theorem 6.1] 
implies that 

q(x, y)~j {g(x + 1, y) + g{x - 1, y) - g(x, y - 1) - g(x, y + 1)) 

is identically zero because it is harmonic with sub-linear growth and attains the 
value 0 at (0, 0). Therefore there exists c\ > 0 such that for all (x, y) € C, we have 

q(x,y) = J {g{x + l,y) + g(x - l,y) - g(x,y - 1) - g(x,y + 1)) + ^Ag(x, y) 

= \ {d(x + 1,2/) +g(x - 1,2/) - 2 g(x,y)) 
x 2 — y 2 

> Cl ( x 2 + y iy 

The first line above follows from the fact that Ag = 0 for all points except (0,0). 

The last line above follows from [Law96, Theorem 1.6.5 (b)]. Therefore itoo(l, 0) > 

2 2 

f E {x,y)ec q( x ’ y) > 52(x,y)eB . jE = oo. Hence s 0 does not stabilize. □ 


We need the following technical lemma for the next example: 


Lemma 6.2. Let a : V —>• R be a configuration in G = (V,E). Assume that 
H : = {x : cr(x) > 1} satisfies \H\ < oo and 

^(u(x) - 1) + < oo. 
x£ H 

Let F C V be such that [F| < oo and 

^(a(x) - 1) + < J^(l - o-(x)) + . 

xeV x&F 


Then a stabilizes. 
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Proof. Let a : H xF —>• [0, oo) be a non-negative function such that YlyeF a(x, y) = 
a(x) — 1 for all x € H and YIxgh a (+ v) — 1 ~ a (y) f° r 1 / € F. The function a 
encodes how to redistribute the excess mass from H to F. 

Let 

g z (x,y) E x (number of vists to y before being killed at z) 

X ' Z V deg (y) deg (y) 

denote the Green’s function normalized with degree. Observe that A f XjZ = 6 Z — S X . 
Therefore the function u := YlxeH y&F a ( x iV)fx,y satisfies a + Au < 1 and u > 0. 
This in turn implies that a stabilizes. □ 

Remark. The condition that \H\, |F| < oo in the above lemma is necessary. The 
following example illustrates this: Consider a probability measure /i on N* = 
{1,2,3,...} and consider the function a y = 1 + 5q — y(x) where y.{x) = /i({x}). 
Then it can be shown that stabilizes if and only if x r( x ) < °°- 

Lemma 6.3. Define C a = {( x ,y) € 1? : x > 0, |y| < ax}. Then the configuration 
s a = rn.lc a stabilizes if \ Tffi < 1 and a € (0,1]. Moreover so(x,y) = ^l{a;>o,y=o} 
stabilizes. 


Proof. The case a = 1 is trivial. Define for a £ (0,1] 


u a ( x ,y) = 


(ax — 


,2\ ^-Ca- 


2(1 +a 2 ) 

To see that so stabilizes, we check that so +Aui < 1. This follows immediately 
from the computation of Aui as 

/ 

1 — a: if x > 0, y = 0 

1 


Am(x,y) = < 


if \y\ < x, x > 0 
if \y\ = x, x > 0 
if x = y = 0 
otherwise. 


A direct computation yields A u a in different regions: An a (0,0) = jrufiF) — T 
If y = 0, x > 0 and all neighbors of ( x,0 ) are in C a (i.e. x > [1/a]), then 


Au a (x,0) = 1- 


2a 


1 + a 2 


rX. 


If y = 0, x > 0 and ( x , ±1) ^ C a , then 

Au a (x, 0) = 


a 2 (l — x 2 ) 


1 + a 2 ' 

If all neighbors of (x, y) are in C a with y / 0, then 

A u a (x,y) = 1. 

If (x, y) € C a with y 0 and one of the neighbors is not in C a , then 


0 < 


1 + a- 


< Au a (x,y) < 1. 



THE DIVISIBLE SANDPILE AT CRITICAL DENSITY 


21 


If (x,y) ^ C a and if all neighbors of (x,y) are not in C a , then 

Au a (x,y) = 0. 

If (x, y) (j C a and one of the neighbors in C a , then at-least 2 of the neighbors are 
not in C a and 

0 < A u a (x,y) < 1. 

Consider v a (x,y) = u\(x + [l/a],y) — mu a (x + [l/a],y). It is easy to check that 
Va> 0ifg^<l. 

If x > 0 we have 

s a (x, 0) + Au a (x,0) = m + (1 - [1/a] - x) - m (1 — a (x + [1/a])) 

/ 2 am \ , _ . 

- 1 + + a 2 ~ 1 J ( x + T 1 / 0 !) 

< l. 

If x > 0, y ^ 0, (x, y) € C a , then 

s a (x,y) + Av a (x,y) = m + l-m = l. 

If V + 0, (x,y) ^ Ca, then 

s a (x,y) + Av a (x,y) < Am<l. 

If x < — [1/a] — 1, then 

s a (x,y) + Av a {x,y) = 0. 

Therefore s a + Av a < 1 for all points except on the finite set {(x, 0) : — [1/a] — 1 < 
x < 0} and s a + Av a = 0 for all x < —1 — [1/a]. Lemma 6.2 implies that s a + Av a 
stabilizes, and therefore s a stabilizes. □ 

We conjecture that the bound in Lemma 6.3 is sharp. 

Conjecture 6.4. Define for a € (0,1] , C a = {(x, y) € 1? : x > 0, \y\ < ax}. Then 
the divisible sandpile s a = mlc a stabilizes if and only if < 1. Furthermore, 
the divisible sandpile so = kx l{ x> o, y=o} stabilizes if and only if k < 1. 

More generally, we have the following problem. 

Open problem 6.5 (Tests for stabilizability). Given s : Tfi —>• M, find series tests 
or other criteria that can distinguish between stabilizing and exploding s. 

7. Finite graphs 

Let G = (V,E) be a finite connected graph with \V\ = n. For a finite con¬ 
nected graph, all harmonic functions are constant: the kernel of A is 1-dimensional 
spanned by the constant function 1. 

Lemma 7.1. Let s : V —> R be a divisible sandpile with s ( x ) = n - Then s 

stabilizes to the all 1 configuration, and the odometer of s is the unique function 
u satisfying s + A u = 1 and min it = 0. 
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Proof. Since A has rank to — 1 and Y1xgv( s ( x ) — 1) = 0, we have s — 1 = Av for 
some v. Letting w = v — mint), we have w > 0 and s + Aw = 1, so s stabilizes. 
Now if to is any function satisfying s + Ato < 1, then 

+ Au)(x) = TO 

xev 

so in fact s + An = 1. This shows that s stabilizes to the all 1 configuration, 
and moreover any two functions to satisfying s + Ato < 1 differ by an additive 
constant. By the least action principle (Proposition 2.5), among these functions 
the odometer is the smallest nonnegative one, so its minimum is 0. □ 

Fix x,z £ V and let f(y ) = 9 ^g{y) be th e function satisfying /(z) = 0 and 
Af = S z — 5 X . (Here g z (x,y ) is the expected number of visits to y by a random 
walk started at x before hitting z). With a slight abuse of notation, we define 
g(x,y) := E 2G y ±g z {x,y). 

Proof of Proposition 1.3. Observe that ExeV 5 ^) = n - Therefore s stabilizes to 
the all 1 configuration by Lemma 7.1, and the odometer u satisfies 

s + A to = 1 

and min u = 0. 

Since A j = 5 Z — 8 X , the function 

v z (y ) := i \ * J^g z (x,y)(s(x) - 1) 

d eg(y) ^ 

has A v z (y) = 1 — s(y) for y / z and 

A v z (z) = y^(s(x) — 1) = 1 — s(z). 

X^Z 

Thus u — v z is harmonic on V and hence is a (random) constant. 

Let to = g ZzeV yZ - Since u — v z is constant for all z, the difference to —to is also 
constant. Recalling that g = g E zS v g z , we have v(y) = gg^y Ex ev 9( x , v)( s ( x )~ 
1). To compute the covariance of the Gaussian vector to, note that 

E[(s(z) - 1 )(s(w) - 1)] = l {z=w} - - 

1 ’ TO 

hence 

E[v(x)v(y)] 

= deg(x)deg(y) - l)(sW - 1)] 

-« ~ k (i>H (l>H)' 

The function K(y) := E«,eV s(«h y) has A K = J2 z , we v F (<** ~ <*«>) = 

0, so K is a constant. The second term on the right is just gg. Letting C 
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d= 1 

d = 2 

d = 3 

d = 4 

d > 5 

E(»7o - ??x ) 2 < V’ d(n,r) : = 

nr 2 

r z log ( 7 ) 

r 

log(l + ?’) 

1 

Emax{77 x : x G Z^} x 

372 - 

n 


log 77 

(log 77) 1 / 2 


Table 1. Statistics of the 


discrete torus Z(( 


n-Laplacian Gaussian field 77 on the 
In the first line, r = ||x|| 2 and the symbol < 
means there is a dimension-dependent constant C d such that E(7/o — 
7/ x ) 2 < C d ip d (ji, r ) for all x G The second line gives the order of 
the expected value of the maximum of the field up to a dimension- 
dependent constant factor. 


be a 7V(0, ^-) random variable independent of v, the Gaussian vectors 77 and 
(v(x) + C) X £V have the same covariance matrix, so 

77 = v + C. 

Since u — v is constant and min u = 0 we conclude that 

u = v — min v = 77 — min 77 . □ 


8. Green function and bi-Laplacian field on 

The rest of the paper is devoted to the proof of Theorem 1.2. Taking Proposi¬ 
tion 1.3 as a starting point, the expected odometer equals the expected maximum 
of the bi-Laplacian Gaussian field 77 , since 

Eu(.t) = E(? 73 ; — min 77 ) = —Emin 77 = E max 77 


where we have used that E ij x = 0. From the covariance matrix for 77 we see that 


nVx-Vy) 2 = Y, (7k 


z,x) 


zev 


\deg(x) 


ff(z,y) V 

deg(y)J 


We will use asymptotics for the Green function g of the discrete torus Z^ to 
estimate the right side. This will enable us to use Talagrand’s majorizing mea¬ 
sure theorem to determine the order of Emax{ 77 x : x € Z^} up to a dimension- 
dependent constant factor. These calculations are carried out below and sum¬ 
marized in Table 1. The table entries give bounds up to a constant factor de¬ 
pending only on the dimension d. For example, the d = 3 column means that 
there is a positive constant C such that E(t 7 o — 77 x ) 2 < C ||x || 2 for all x e Z^ and 
G _ 1 n 1 / 2 < E max 77 <r,d 


Remark 1. For the rest of our work, we identify the discrete torus Z^ with 
(Zn (—ri/2, n/2]) d which in turn is viewed as a subset of M rf . For x G Z^ and 
1 < p < 00 , we denote by ||x|| p the p-norrn under the above identification. Note 

that for standard graph distance do on Z^, we have dc( 0 5 x ) = || x lli- 
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8.1. Fourier analysis on the discrete torus. In this section, we derive a for¬ 
mula for K(r/o — i?x) 2 - We begin by recalling some basic facts about Fourier analysis 
on the discrete torus and the spectral theory of the Laplacian. We equip the torus 
Zj( with normalized Haar measure g, (in other words the uniform probability mea¬ 
sure). Consider the Hilbert space H = L 2 (Z^ ;/ u) of complex valued functions on 
torus with inner product 

if, 9) = f fgdg = \ Y f( x )9(x)- 

ht f-' 


We identify the Pontryagin dual group Z(( with Z^ as follows. For any a. £ Z d , 
the map x 1 — > exp (i27rx • a/n) gives the corresponding element in Z(( (The dot 
product is the usual Euclidean dot product in M”). We denote this character by 
% a . Recall that {x a : a £ Z^} forms an orthonormal basis for 7~L. Moreover each 
Xa is an eigenfunction for the Laplacian A with eigenvalue 

d 



i =1 


Thus the Laplacian A : 77 —>• 77 is a non-positive, bounded operator. Moreover 
A a = 0 if and only if a = 0. Laplacian A is a self-adjoint operator, that is 


(/ 1 ,A/ 2 ) = (A/ 1 ,/ 2 ) (17) 

for all /i ,/ 2 € 77 (See Remark 2). We denote by g x ( y) = g{ y,x). Recall that 

A g x = 2d (~dXo ~ . (18) 

Denote by g^( a), the Fourier coefficient (<? x , xa}- Since the function x i-a g x (y) 
is harmonic, it is constant. This implies that there exists L > 0 such that 

ftc(O) = n~ d Y 5 x(y) = L (19) 

y&A 

for all x£Zj . For a / 0, we have 

■^a <?x( a ) = A a (5x;X a) = ( 5 x 1 AXa) = (^9x,X a) = — 2d(5 x , Xa) = ~2dn X—a(x). 

( 20 ) 

For the above equation, we used Ay a = A a x a , equations (17), (18) and (xo- Xa) = 
0. By Parseval’s theorem and equations (19), (20), 


E(t ? 0 - ??x ) 2 


i 2d ) 2 Y (s^O) -£f(z,x)) 2 

(2d)~ 2 n d (g 0 - g x ,g 0 - g x ) 

(2 d)~ 2 n d Y l5o( z ) -fe ( z )| 2 
zez^ 

\ F nA*) 


(21) 
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where 


F n ,d(x) ■= n d ^2 


sm 




«3S\{0> (E?.,sin 2 (=?))' 
Remark 2. In M c / we have the Green’s second identity 


( 22 ) 


J /iA / 2 dx = - J 


for all /i, / 2 € 


fiAf 2 dx = - / V/ 1 .V /2 dx = / /iA / 2 dx 


. Similarly, in our discrete setting we have 


(/1.A/2) = X] (-M x ) “ /2(y))(/2(*) - /2(y))fc(x,y) = (A/1,/2) 


x,yez^ 


where /c(x, y) = and x ~ y if x and y are neighbors in Z^. See for instance, 
[Sal96, Lemma 2.1.2] or [LPW09, Lemma 13.11] for a proof. 


Our task now is to estimate the expression -F n ,d(x). Henceforth we assume 
that x/0. To study the quantity max x>yeZ d (r/ x — %) as n goes to 00 , we want 
to estimate E(r/o — ?? x ) 2 with d fixed and n large for different values of x. We 
approximate -F n ,d( x ) by an integral of a function over W l . For w £ and r > 0, 
we denote by B oa ( w, r) the open ball with center w and radius r under supremum 
norm, that is 

Hoo(w, r) = {y € R d : ||y - w]]^ < r}. 

We denote the indicator function of the ball B^z/n, l/(2n)) by / z>n : R d —>• (0,1}, 
that is I x>n = l Boo ( z / n ,i/(2 „))• Define the function G n)d)X : -A- R 


G 


n,d,x — 


E 


Sill 


(™) 


fet 1 S™ 2 ( 


2 /z,n ■ 


Since the cubes B 00 (z/n, 1 /( 2 n)) are disjoint with volume n d , we have 


Fn,d 



G n ,d, x ( y)dy. 


By triangle inequality, we have 


(23) 


(1 + yfd) 1 ||z/n || 2 < ||y || 2 < (1 + Vd) ||z/n || 2 (24) 

for all z € Z^\ { 0 } and for all y € B^z/n, 1 /( 2 re)) under the usual identification 
from Remark 1. We will estimate the function G nj d,x using the function H n ^,^ ■ 
—>• R defined by 

■ 2 / TTXZ \ 

Hn,dA y) = E ,| ,| 4 n ^ n(y)• 

zezd\{ 0 } l|y|| 2 


More precisely, we have the following lemma. 
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Lemma 8.1. Fix d G N*. There exist positive reals c\,C\ such that 

ci#n,d,x(y) < Gn,d, x (y) < CiH n4 ^{ y) (25) 

for all n G N* ; for all x G Z^ \ {0} and for all y G R d . 

Remark 3. We will use Ci for large constants and c* for small constants. Here and 
in what follows, all constants are allowed to depend on d but not on x e Zf \ {0} 
or n. 


Proof. The idea is to use the estimate 


— \t\ < |sint| < \t\ 

IT 

for all t G [ —7 t/2 , tt/ 2] . Thus there exists a constant C 2 > 0 such that 

2 




|z/n ||2 < ^y^sin 2 ^< C 2 ||z/n| 


(26) 


for all z G Z^ and for all n G N*. By (24), we have 

(i + ^-4 < (1 + Vd)*^y\ 

ll z /™|| 2 IMI2 ll z /«ll 2 

for all z G Z^\ {0}, for all y G and for all n G N*. Combining equations 
and (27) gives (25). 


(27) 

(26) 

□ 


By (21), (23) along with integration of (25) over the variable y, there exists 
ci , C\ >0 such that 


cid 2 / 77 n ,d,x(y) dy < E(770 - r/ x ) 2 < Cid 2 / x (y) dy 


(28) 


for all n G N* and for all xGZj \ {0}. 

By (28), it suffices to estimate f Rd H n ,d,,x(y) dy. Observe that the support of 
H n d.x satisfies 


Support(77 n ^ x ) C B 2 { 0 , Vd) \ B 2 { 0 , 1/(2 n)) 


for all d, n 

G N* and for all x G Z^ 

where H 2 

denotes open ball with respect to 

Euclidean 

norm in R d . 




8.2. Upper bounds. Define ipd by 






9 

nr 

if d = 1 




r 2 log (y) 

if d = 2 



ifd{n,r) := < 

r 

if d = 3 

(29) 



log(l + r) 

if d = 4 




1 

if d > 5. 



for all n G N* and all r > 0 along with ^(n, 0) := 0 for all d,n G N*. The upper 
bounds for E(??o — i?x) 2 is summarized in the following Proposition. 
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Proposition 8.2. For each d E N*, there exists C d > 0 such that 

E(r? 0 - i/x) 2 < C d if d (n, ||x|| 2 ) (30) 

for all n E N* and for all x E ; where if d is defined by (29). 

Proof. By (28) it suffices to find upper bounds for f Rd H njdjX (y) dy. The strategy 
to establish upper bounds for f Rd H n ,d,x(y) dy is to split it into two integrals 

aS h* = /l/(2n)<||y|| 2 <Vd/||x|| 2 +/v / d/||x|| 2 <||y|| 2 <Vd- N ° te that b ° th tlle inte S rals 
are over non-empty annuli since l/(2n) < \fd/{A ||x|| 2 ) < \/d/||x|| 2 < Vd for 
all x E \ {0} with the identification from Remark 1. Using Cauchy-Schwarz 
inequality and the bound |sint| < t, we have |sin(7rx.z/n)| < 7r || x || 2 ||z/n|| 2 . This 
bound competes with the trivial bound |sin(7rx.z/n)| <1. It will become clear 
that up to constants, the first bound is better for the first term and the trivial 
bound |sin(7rx.z/n)| < 1 is better for the second term. 

For the first integral we use the bound |sin f| < |f| and Cauchy-Schwarz inequality 
to obtain 

H n ,dA y) < 

zez£\{o} 

By (24), we have ||z/n|| 2 I x , n (y) < (1 + Vd) 2 ||y || 2 / z ,„(y). Therefore, we obtain 


2 II ii2 
7T X 


2 l|Z/n| -/ Zi „{y). 


|y|l2 


(1 + Vd) 2 ir 2 


I II2 

| y || 2 


-4,n(y)- 


H n ,d,V y) < 

zGZ^\{0} 

for all n, d E N*, for all y E and for all xEZj\ {0}. Hence, we have 


(31) 


/, : = 


-l 


1 /(2n)<||y|| 2 <\/d/||x|| 2 
2tT d ' 2 


H n ,d, x ( y) dy < (1 + Vdfir 2 


\x\\ 2 Ud-l 


L 


\di/l|x|i2 r < 
1/(2 n) 


d-1 


dr 


(32) 

where u d _i = Y[d/ 2 ) — l)-dimensional surface measure of unit sphere § rf_1 


m 


For the second integral we use the bound |sini| < 1, to obtain 


h : = 


Hn,dVy)dy < 1 


rVd 

Vd/ |b 


„d -1 


■ dr. 


(33) 


d— 5 


' v/rf/||x|| 2 <||y|| 2 <v / rf 

Combining equations (32) and (33), we obtain 

r f'/d/ ||x|| 2 

/ Hn,dAy)dy < (l+\/d) 2 7T- ||x|||/ r d_3 dr+w rf _i 

J l/(2n) 

(34) 

The desired upper bounds on E(r/o — i? x ) 2 for all dimensions follow from (34) along 
with (28). □ 


rVd 

'Vd/ ||j 


dr. 


Remark 4. The terms /] and I 2 correspond to the energy (square of 2-norm) 
of the low and high frequency oscillations of the function go — g x respectively. 
For d = 1,2, the term I\ dominates 1 2 - For d = 3, both I\ and I 2 are of the 
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same order. For d > 4, the term I 2 dominates I\. Hence our approach to obtain 
matching lower bounds in the next subsection is as follows: For d = 1,2, we obtain 
lower bounds on lower frequency terms and for d > 3 we obtain lower bounds on 
higher frequency terms. 

It is well know that the Gaussian field r) induces a Hilbert space on Z(( given 
by the distance metric 

<^(x,y) := (E(r7 x -?7 y ) 2 ) 1/2 . 

The upper bounds on d v provided by Proposition 8.2 transfers to upper bounds 
on IEsup xgZ d ?? x . The main tool to transfer bounds is Dudley’s bound [Tal05, 
Proposition 1.2.1] described below. There exists L > 0 such that 

OO 

E sup < L E 2 k/2 e k (35) 

k=0 

where e k = inf sup tgZ d d v ( t, T k ) and the infimum is taken over all subsets T k C 
with \T k \ < 2 2k . 

Proposition 8.3. For each d € N*, there exists C d > 0 such that 


E sup r/ x < C d fi d {n) (36) 

xez£ 

for all n € N* , where <f> d is defined by (2). 

Proof. Let do denote the standard graph distance on the torus ZBy choosing 
a submesh of appropriate cardinality the following statement is clear: For any 
d £ N*, there exist C d) \ > 0 such that for any n > 2 and for any 2 < m < n d , 
there exists a set S m C TLi with |5 m | = m such that 

Ti 

sup d G {t, S m ) = sup inf d G (t,s) < C d ,i—r-n- (37) 

tez^ tez ^m / 


For each d € N* by Dudley’s bound (35), (37) and Proposition 8.2, there exists 
C d) 2 , C d) 3 > 0 


[log log n d \ 

E max ri x <C d 2 2 k ^ 2 

xez " ’ to 


'fid n, 


C d ,in \ 
2 ( 2 k l d ) ) 


I 1/2 


< C d , 3 fi d {n) (38) 


The second inequality above follows from a straightforward case by case calcula¬ 
tion. □ 


8.3. Lower bounds. Next, we prove matching lower bounds on E(?]o — ??x) 2 - For 
dimensions d = 1,2, we estimate F U)d directly. 
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d = 1: We use the bound |t| > |sin£| > ^-\t\ for all \t\ < 7 r /2 to obtain 


sin 2 ZE 2 S 

F n , i(x) > rT 1 . 4 n > 47 r _4 ?r ||x||| 
sin - 


for all n G N* and for all x G Z n \ {0}. Hence there exists ci > 0 such that 

IE(r/o - f?x) 2 > cin ||x||| (39) 

for all n € N* and for all x G Z n \ {0}. 

d = 2: Let S'*. C Z 2 denote the sphere with center 0 and radius k in the supremum 
norm, that is S*, = {y G Z 2 : ||y= k}. For x G M 2 , we define H x = {y G Z 2 : 
|x • y| > ||x|| 2 ||y|| 2 /V2}. It is easy to check that \Sk\ = 4 k and \Sk (~l Lf x | > 2k for 
all k G N* and for all xGl 2 . Let a G (1, \/2). If ||z || 2 < by Cauchy-Schwarz 

inequality we have |x • z/n| < a -1 . We need the inequality n\t\ > |sin tt£| > f3\t\ for 
all |f| < a -1 where /? = a sin( 7 ro; _1 ). Putting together the above pieces, we obtain 


F n , 2 (x) > 


J2 E^ 27r " 4 ii z ii2 4 i x 

k= 1 zeSfc 

„ 2 [«%H"/ 4 J 

^4 H X ll2 Z. H Z 

k =i z,eS k nH x 

„ 2 [n&hwu 

£iiwl E E r 

fc=l zeS^D-Hx 

„ 2 [^pl-/ 4 ) 

E *- 


> ^ II ||2 

a pH. 


. /3 2 „ i 2 

- 2G l|X " 2 


. p 1 . ||2 , 

a 2plHI 2 i»s 
„ /? 2 | l|2 . 


A |n/4j ) + 1 


n n 

n n A- 
a ||x || 2 4 


Since ||x || 2 < n/y /2 for all x G Z 2 , we have the desired lower bound by using 
a < y/2. That is, there exists C 2 > 0 such that 


E(??o - i/x ) 2 > c 2 11 x 112 log f 


for all n G N* with n > 4 and for all x G Z 2 \ {0}. 
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d > 3: For dimensions d > 3, we will approximate H n ,d, x by its almost everywhere 
point-wise limit iioo,d,x defined by 


-^oo,a!,x(y) 


sin 2 (7rx • y). 


L Rc 


do,i/2) (y)- 


for y / 0 and 0 otherwise. Therefore we would like to estimate integrals of the 
form 


/ 

J ri<|| 


sin 2 (7rx • y) 


,<r 2 


dy = 


rr 

J r i 


A— 5 


Sd(r 


X 


i dr 


(41) 


where Sd(t) = fgd-i sin 2 (7rtyi) Vd-i{dy) and Vd-i denotes the surface measure in 
S d—1 . We will need the following lower bound for Sd- 


Lemma 8.4. Fix d E N* with d > 3. Then for all e > 0, there exists 5 > 0 such 
that Sd(t ) > 5 for all t > e. 

Proof. Since Sd : M —>• M is a continuous function with Sd(t ) > 0 for all t 0, it 
suffices to show that liminf Sdft) > 0. By [BGMN05, Corollary 4] (See Remark 

t—> OO 

5(b)), 

r l fl/2 

Sd(t) =Cd (1 — x 2 )( rf - 3 )/ 2 sin 2 (7rtu;) dw > Cd2^~ d ^ 2 / sin 2 (7rtu.!) dw 

J- 1 J- 1/2 

where is a constant that depends on d. Since lim^oo sin 2 (7rtu;) dw = 1/2, 
the conclusion follows. □ 


Remark 5. (a) Recall that we used the point-wise bound |sint| < 1 to obtain 

upper bounds on I 2 . We want to somehow reverse that inequality to obtain 
corresponding lower bounds. Although the reverse inequality |sint| > 5 is 
not true for any 6 > 0 in a pointwise sense, it is true in an average sense. 
That is the content of Lemma 8.4. 

(b) [BGMN05, Corollary 4] implies the following striking result in geomet¬ 

ric probability: Let d > 3. For a uniformly distributed random vector 
y = ( 2 / 1 , 2 / 2 , • • • , Ud) hi the (d — l)-dimensional unit sphere in W l , 

the projection (yi, 7 / 2 ,..., yd- 2 ) is uniformly distributed in the (d — 2)- 
dimensional unit ball M d ~ 2 = LL(0,1) in M rf ~ 2 . 

(c) Since lim^^oo H n = L^oo,d,x almost everywhere, one might wonder if 
we can prove matching lower bounds for I 2 using dominated convergence 
theorem. This approach gives a lower bound as n goes to 00 but with both 
d and x fixed. However we want lower bounds with fixed d and with both 
n and x varying. Hence there is a need to quantify this convergence as 
both n and x varies. We fulfill this need in Lemma 8.5. 


One can easily check that lim^^oo i4 n ^ x = H^ d,^ almost everywhere. We need 
the following quantitative version of this convergence. 
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Lemma 8.5. Fix d £ N*. For any e > 0, there exists positive reals S,N such that 

sin 2 ( 7 rx • y) 


Hn,d,x. (y ) 


I l|4 

|y || 2 


< 


1 n -i 

|y || 2 


for all n > N, for all x € Z(( \ { 0 } with 
B 2 (0, 1/4) \S 2 (0,1/(8 ||x|| 2 )). 


|x || 2 < Sn and for almost every y € 


Proof. Note that, we have the inclusion 


B^O, 1/(2n)) C S 2 (0, V(8 ||x|| 2 )) C B 2 { 0,1/4) 

for all n > Ay/d and for all x £ h d \{0} with || x|| 2 < n/Ay/d. We use ||x|| 2 > 1 and 
the comparison of norms ||w|| 2 < x/dllwH^ to prove the above inclusions. The 
function y i-A sin 2 (7rx • y) has gradient bounded uniformly in 2-norm by 7r ||x|| 2 . 
Hence we have 


Hn,d,x (y ) 


sin 2 (-/rx • y) 


< ^ IMI 2 


Vd 

n 


y 


-4 

2 


for all n > Ay/d and for all x £ Z^ \ {0} with || x 11 2 n/Ayfd, and for almost every 
y £ i? 2 ( 0 ,1/4) \ B 2 (0, 1/(8 ||x|| 2 )). The choice 6 = min ^ 3 ) and N = Ay/d 
satisfies all the requirements. □ 


We put together the above pieces to obtain the following lower bound for d > 3. 
Lemma 8 . 6 . Fix d> 3. There exists positive reals 5,N,Cd such that 

r 


H n ,d,Ay) dy > c d 


r d 5 dr 


'1/(8||x|| 2 ) 


for all n > N and for all x £ Z^ \ {0} with ||x|| 2 < Sn. 

Proof. By Lemma 8.4, there exists e\ > 0 such that Sd(t) > 2ei for all t > 1/8. 
By Lemma 8.5, there exists positive reals 6, N such that 


H n ,d, x ( y) 


sin 2 (7rx • y) 


I l|4 

|y || 2 


< 


£1 

Ud -1 lly || 2 


for all n > N, for all x £ Z^ \ {0} with ||x || 2 < Sn and for almost every y £ 
£ 2 ( 0 ,1/4) \B 2 (0,1/(8 ||x|| 2 )). 
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Combining the above observations, we have for all n > N and for all x E Z^\{0} 
with ||x || 2 < fin 


H n ,dAy) dy> 

J i 

>-l 


l/(8||x)|| a <||y|| a <l/4 


Hn,d,*(y) dy 


(Hoo,d, x (y) ~ e ^d\ lly || 2 4 ) dy 

V(8||x)|| 2 <||y|| 2 <l/4 
/•1/4 

= / r d_ 5 (s d (r ||x|| 2 ) -ei)dr 

•' 1 /( 8 II X II 2 ) 

/*!/4 

> ei / /" 5 dr. □ 

*' 1 /(8|I x II 2 ) 

We now establish the following lower bounds corresponding to the upper bounds 
in Proposition 8.2. 


Proposition 8.7. For each d E N*, there exists positive reals 5d,Nd,Cd such that 

E(r?o - Vx ) 2 > c d ^d(n, ||x|| 2 ) 

for all n > Nd and for all x£Zj\ { 0 } with ||x|| 2 < fidn, where ifd is defined by 
(29). 

Proof. The cases d = 1,2 follow from (39) and (40) respectively. The case d> 3 
follows from Lemma 8.6 along with (28). □ 

Remark 6 . The condition ||x|| 2 < fidn that appears in the lower bound for the case 
d > 3 is somewhat unsatisfactory. We believe that the lower bound is true without 
any such an additional condition. However the lower bounds in the present form 
are good enough for our main application. 

Next, we obtain lower bounds matching the upper bounds in Proposition 8.3. 
We start by recalling notation and setup for Talagrand’s majorizing measure 
[Tal05, Theorem 2.1.1]. We consider centered multivariate Gaussian random vari¬ 
ables ( rjt)teT indexed by a set T with cardinality |T|. An admissible sequence {^4^} 
is an increasing sequence of partitions of T such that |^| < 2^ . Here “increasing 
sequence” refers to the fact that every set in A n+ \ is contained in a set in A n . 
We denote by Ak(t) the unique element of A n that contains t E T. Recall that 

d v (ti,t 2 ) = (E(? 7 t i — rfe ) 2 ) 1 / 2 denotes the Hilbert space metric induced by ( r]t)teT■ 
We define the function 

OO 

72 {T,drf) = inf sup V'diam r? (A n (t)) 
teT k =o 

where diam, ; denotes the diameter in the d v metric and the infimum is taken over 
all admissible sequences. The majorizing measure theorem [Tal05, Theorem 2.1.1] 
states that there is some universal constant L for which 

y 7 2 (T, d) < Esup rj t < Ly 2 {T,d). 
teT 


( 42 ) 
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Proposition 8.8. For each d £ N*, there exists Cd > 0 such that 

E sup ?? x > Cd</>d(n) (43) 

for all n E N* ; where fid is defined by (2). 

Proof. The strategy is to use the lower bound for d v given by Proposition 8.7 along 
with (42) where T = {x £ Zj : ||x || 2 < hdn} and Esup igZ d rjt > Esup igT i?t. Note 
that it suffices to show (43) for large enough n, i.e. n> Nd f° r some fixed Nd- 
For d = 1,2,3, by (42) we have 

E sup r]t > L~ 1 inf sup diam^ (Ao (t) ). (44) 

tezi teT 

Since |Ao| < 2 , we have sup tgT diamG(Ao(t)) > con for some Co > 0. Therefore by 
Proposition 8.7 along with (44) we obtain the desired result. 

For d = 4, by (42) we have 

E sup rjt > L inf sup 2 k / 2 diam^A^f)). (45) 

tez£ teT 

where k = [log 2 log 2 |F|J —1. This gives 2 fc / 2 > co^/\ogn for some co > 0. Moreover, 
k = |_log 2 log 2 |T|J — 1 and \Af\ < 2~ implies that at least one of the sets A^if) has 
cardinality greater than or equal to which in turn implies diamG>(j4fc(t)) > 

c\s/n for some ci > 0. By Proposition 8.7, we obtain diamG( J 4fe(t)) > C 2 \fiogn 
for some c 2 > 0 and for large enough n. The conclusion for d = 4 then follows 
from (45). 

The case d > 5 is a direct consequence of Sudakov minoration ([Tal05, Lemma 

2 . 1 . 2 ]). □ 

Proof of Proposition 1.2. The upper and lower bounds follow from Propositions 
8.3 and 8 . 8 . □ 
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